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Abstract 

Within the framework of the Brueckner theory, the off-shell behaviors of the mass operator 
M{k,uj) = V{k,uj) + iW{k,u}), i.e., its dependence upon the momentum k and upon the nucleon 
frequency co, are investigated by including nuclear three-body force (TBF). The first two terms 
of the hole- line expansion of the mass operator are taken into account. The TBF effects on their 
off-shell properties are discussed. A comparison is made between the on-shell and off-shell values 
of Ml. The nucleon spectral function and nucleon momentum distribution are also calculated, and 
the calculation shows that they are hardly affected by the TBF effect at the saturation density. At 
a high density two times greater than the saturation density, inclusion of the TBF may lead to a 
visible effect on the spectral function and may enhance the depletion of the hole states. 
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I. INTRODUCTION 



The probability of removing a particle with momentum k from a target nuclear system, 
leaving the final system with excitation energy u, is reflected by the nucleon spectral function 
S{k,uj). In the free Fermi gas model, the spectral function can be written as S{k,u}) = 
6{u) — h^k'^/2'm). However, many-body correlations among the nucleons, induced by the 
nucleon-nucleon {NN) interactions, broaden the peaks of the Fermi-gas spectral function 
and decrease their strengths. It has been shown that, most of this decrease of the strength 
of the single-particle (s.p.) states with respect to the standard mean field estimates |l[ is 
due to the A^A^ correlations {2, 3|, whose effects can be most accurately studied in infinite 
nuclear matter. Therefore, a microscopic calculation of the nucleon spectral function in 
nuclear matter is of special interest since it may play an important role in understanding 



QQ. 



the nature of the NN correlations, especially the short-range and tensor correlations 

The interest in the spectral function has also been raised by the treatment of the off-shell 
effect in transport theory, which provides a generalized theoretical framework to describe 
the time evolution of heavy ion reactions. The use of the quasiparticle approximation (QPA) 
in transport theory puts the nucleon on the mass shell, neglecting not only the finite decay 
width of the particles, but also the width of the nucleon spectral function. However, the 
on-shell quasiparticle limit should not be adequate for particles with short lifetimes and/or 
high collision rates as recognized previously 6|,Q]. Therefore, there have been attempts to go 



beyond the QPA, such as transport formulations for quasipartic 
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es with dynamical spectral 



2l|. The nucleon spectral 



functions |8l-ll0|. and extend the extensively applied models 
function plays an important role in the implementation of the off-shell effects in a transport 
theoretical treatment of heavy-ion and other nuclear collisions. 

Experimentally, the information about the nuclear spectral function and the effect of 
A^A^ correlations in nuclear systems can be extracted from the (e, e'p) and proton-induced 



knockout reactions 



flfl 



22l-l26|. Theoretically, the nuclear short-range correlations and 



the spectral function in nuclear matter have been investigated extensively by using vari 



ous microscopic nuclear many-body approaches, such as the Green function theory 27H34|. 



the correlated basis function method 
framework 



35 




the extended Brueckner-Hartree-Fock (EBHF) 



and the in-medium T-matrix approach 



readers to Refs. [4, . Within the framework of the Brueckner theory, the nucleon spectral 



44j-|46l|. For a review, we refer 
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function in symmetric nuclear matter has been studied in Ref. 4l| by adopting a finite-rank 
representation of the reahstic Argonne VIA NN interaction, without taking into account 
any three-body force (TBF) effect. In Refs. SJ, |43|, the neutron and proton spectral func- 
tions in isospin-asymmetric nuclear matter have been explored using the BHF approach 
and the Green function theory, respectively. It is well known that inclusion of TBF in the 
nonrelativistic Brueckner theory is crucial for reproducing the nuclear saturation properties 
and for better describing the s.p. properties, such as the momentum dependence of the 



nucleon s.p. potential 



47 



50|. Recently, the TBF effect on the spectral function in nuclear 



matter has been investigated within the framework of the in-medium T-matrix method in 
Ref. H, where the TBF adopted is the Urbana TBF [5l|. In that paper, the authors have 
shown that the TBF effect on the spectral functions is quite small at low densities around 
and below the saturation density and that noticeable modification of the spectral functions 
is realized only for high densities well above the saturation density. One of our purposes in 
the present paper is to investigate the possible impact of a microscopic TBF on the nucleon 
spectral function within the framework of the extended BHF approach. 

The spectral function is closely related to the mass operator M{k,u)), whose off-shell 
behavior is also our concern in the present paper. The off-shell mass operator plays an 
important role in the dispersion relation to the nuclear mean field 52|] and in the discussion 



of y-scaling in inclusive electron scattering 



53|-|55|. In Ref. 



41| . the properties of the off- 



shell mass operator have been obtained within the Brueckner theory in the absence of any 
TBF. Therefore, the other purpose of the present paper is to reveal the TBF effect on the 
off-shell mass operator discussed in Ref. 41|. 

The paper is organized as follows. In Sec. II, we give a brief review of the mass operator 
and spectral function, i.e., their definitions and physical interpretations, as presented in Ref. 



4l| . We also provide a simple introduction to the Brueckner theory and the microscopic 



TBF adopted in our calculation. In Sees. Ill and IV, we focus on the real and imaginary 
parts of the off-shell mass operator. We study the fc-dependence of M{k,u) for the two 
typical energies, u = 20 MeV and u = 160 MeV, at 0.34 fm~^. We also calculate the 
cj-dependence of M{k,u) for k = jkp and k = jkp- The off-shell results of Mi{k,u) are 
compared with the on-shell ones. The TBF effect on the k and w-dependence of M{k, uj) is 
discussed. In Sec. V, we calculate the spectral function and investigate the TBF effect on 
its w- dependence. In Sec. VI, a summary is given. 
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II. FORMALISM 



A. The mass operator and the spectral function 

The Green function in the energy-momentum representation is given by G{k,uj) = [u — 
k^/2m - M{k, uj)]~\ where M{k, oo) = V{k, oo) + iW{k, u) is the mass operator that can be 



identified with the mean field felt by a nucleon in a nuclear system, 
parts of the mass operator are connected by the dispersion relation 



he real and imaginary 



4l|: 



V{k, oo) = lim V{k, uj) + - r El^li^ doo' . (1) 

i^^oo TT oo' — OO 

The spectral function is given by 

^ ' ' n\[uj-ky2m-V{k,uj)]^ + [W{k,uj)]y ' 
and it fulfills the sum rule 

/oo 
S{k,u)du = l. (3) 
-oo 

The occupation probability n{k) is related to the spectral function by 



n{k) = I S{k,oo)du (4) 

and 

n{k) = 1-1 S{k,u)doo . (5) 



The Fermi energy up fulfills oop = kp/2m + V{kF, oop)- For a system of A nucleons, S{k, E*) 
measures the probability density of finding the residual {A — l)-nucleon system with exci- 
tation energy E* = up — u}{oo < oop) after removing a nucleon with momentum k from the 
ground state, or the probability density of finding the residual {A + l)-nucleon system with 
the excitation energy E* = oo — oop{oo > oop) after one has added a nucleon with momentum 
k to the ground state. 



B. Brueckner theory with a microscopic TBF 

The starting point of Brueckner calculation of nuclear matter properties is to obtain the 
Brueckner reaction matrix G{oo) by solving the Bethe-Goldstone (BG) equation 

Gioo) = v.. + v.. "^'f ^(--) ' 
^^uj - t[ki) - t[k2) + iri 
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where ki and ^2 are momenta of the two involved nucleons. Q{ki, ^2) = [1 """-(^i)]!! — ''^(^2)] 

is the Pauh operator which prevents two intermediate nucleons from being scattered into 

occupied states, u is the starting energy. The single-particle energy e{k) satisfies the on-shell 

relation e{k) = A;^/2m + f/BHF(^); where the auxiliary potential U-BRpik) is the single-particle 

potential at the BHF level and it is defined as UBUpik) = '^R.e{kk'\G{e{k) + e{k'))\kk') a. 

k' 

The subscript A denotes antisymmetrization of the matrix element. The continuous choice 



other than gap choice is adopted when solving the BG equation to obtain the G-matrix 



56| 



Extension of the Brueckner-Bethe-Goldstone (BBG) theory to include TBFs can be found 



in Refs. 
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48l | . In this paper, we choose the microscopic TBF which is based on the meson- 



exchange current model proposed by P. Grange et al. [57| and reduced to an equivalent 
effective two-body force V^^ via an average with respect to the third-nucleon degree of 
freedom. The effective force V^^ in r space reads 

n 

X [1 - ^('^23)]^3(^iy2>^3l^l>^2y3) 

X0„(r-'3)[l-^7(ri3)][l-r7(r23)] , (7) 

where the wave function (pn denotes the single nucleon wave function in free space. The 
realistic A^A^ interaction Vnn in the BG equation is the sum of the Argonne Vis (^^18) 
two-body interaction and the effective two-body force ^3*^^, as described in Refs. 47|, |48 |. 



Since ri{r) in expression (j7]) is the so-called defect function j57|, l58|| corresponding to the 
G-matrix, V^^ should be recalculated along with the G-matrix in each iteration of our BHF 
procedure to ensure self-consistency of the BG equation. 

In the spirit of Brueckner theory, the first two terms of the hole-line expansion of the 
mass operator are the BHF approximation Mi{k, u) and the Pauli rearrangement correction 
M2{k,uj). They are represented by the diagrams of Fig.l, and their expressions read : 

Mi(fc,u;) = ^ (A;/i|G[t^ + e(/i)]|A;/i)^, (8) 

h<kp 

1 >^ \{lm\G[e{l)+e{m)]\kn)A\' . . 

l,m<kF,n>kF ^ ' ^ ' ^ ' 

Their off-shell values can be calculated as long as the G-matrix is obtained. 
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FIG. 1: Representation of the first two terms of the hole-line expansion of the mass operator. The 
thin lines represent either particle (upward-pointing arrows) or hole (downward-pointing arrows) 
momentum states. The thick lines show the values of the nucleon momentum k and frequency w. 



III. FREQUENCY DEPENDENCE OF THE MASS OPERATOR AT FIXED MO- 
MENTUM 

A. Real part of the off-shell mass operator 

The calculated cj- dependence of V\{k,uj) = ReMi{k,uj) and V2{k,uj) = ReM2{k,u) is 
shown in Fig. [2]for the two densities of p = 0.17 fm~^ and p = 0.34 fm"^, respectively. Two 
fixed momenta (fc = ^kp and k = ^kp) are selected. The quantity ep is the calculated value 
of the single-particle energy e{k) at the Fermi momentum: ep = e{kp) = /c|,/2m + f/BHF(^F)- 
As one can see in Fig. [2], the quantity Vi{k,uj) is attractive for u < ep and its attraction 
increases as a function of frequency u in the region of u < ep, while V2{k,u) is repulsive 
for uj > ep and its repulsion decreases with increasing u in the region of u > ep. The TBF 
effect on their cu-dependence is also reported in this figure. Inclusion of the TBF in our 
calculations hardly affects the dependence of V2{k,co), but tends to reduce the attraction 
of Vi{k,u) well below ep and enhance its attraction as u is much larger than ep. At the 
saturation density of 0.17 fm~^, the TBF effect on Vi{k,u) is weak enough to be neglected 
in the vicinity of ep. However, the TBF effect gets much stronger at high densities. As a 
result, the TBF-induced reduction of the attraction of Vi{k,u) well below ep is obviously 
seen at two times the saturation density 0.34 fm~'^, as revealed in the right panel of Fig. [2l 
At high densities, the TBF effect on V2 turns out to be rather small. At p = 0.34fm~^ and 
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k = 2.1fm~^, inclusion of the TBF may enhance shghtly the repulsion of V2. 

Besides, it is worth noticing that the distinct deviation of the curve with open squares 
from that with filled squares when u—ep is above 150 MeV. The deviation appears regardless 
of the density value, which indicates that one should account for the TBF effect carefully in 
the high-energy domain. 



B. Imaginary part of the ofF-shell mass operator 

Figure E] shows the dependence of Wi{k,ijj) and W2{k,ijj) upon the difference u — ep- 
One important feature of the two components is that Wi{k,uj) vanishes for u < ep and 
W2{k,uj) vanishes for uo > ep- Moreover, W2{k,u) also vanishes for large negative u. At the 
saturation density of 0.17 fm~^, the calculated Wi{k,u) including the TBF contribution is 
very close to its values without including the TBF contribution in the energy domain ranging 
from ep to approximately 300 MeV. However, as density increases to 0.34 fm~^ where the 
TBF effect becomes strong, inclusion of the TBF leads to a faster increase of the attraction 
of Wi{k,u)) with increasing frequency u as compared to the result without the TBF effect. 
At p = 0.34fm~^ and k = 2.1fm~^, inclusion of the TBF may lead to a sizable enhancement 
of the attraction of W2. 



C. Comparison with on-shell values 

In Fig. m we compare the off-shell values of Wi{k,u) and W2{k,uj) with their on-shell 
values. Although our calculations are done at a higher density of 0.34 fm~^ and in the 
presence of the TBF, the results plotted in Fig. H] are similar qualitatively to those in Fig. 12 
of Ref. 4l|, regardless of the magnitude. Therefore, the analysis and conclusion in Ref. 
4l| remain valid. That is to say, on the one hand, Wi{k,uj) and W2{k,uj) are symmetric 
with each other only in the vicinity of the Fermi energy; on the other hand, the assumption 
in the simplest version of the dispersion relation approach for the nuclear mean field [ i.e., 
the w-dependence of Wi{k,u) is approximated by the e-dependence of the on-shell Vri(e)] 
is only justified qualitatively. 
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FIG. 2: (Color online) Dependence of Vi{k,uj) and V2{k,uj) upon to — ep for the two densities of 
p = 0.17 fm~^ and p = 0.34 fm~'^, and for the two fixed momenta of A; = jkp and k = jkp. The 
curves with open squares and open triangles have taken into account the TBF contribution. 
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FIG. 3: (Color online) Dependence of Wi(k,uj) and W2{k,uj) upon oj — ep for the two densities of 
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FIG. 4: (Color online) Comparison between the e-dependence of the on-shell M^i(e) and W2{e) 
(solid curve) with the w-dependence of Wi{k,uj) and of W2{k,iLi), for k = 1.3 fm^-*^ (filled squares) 
and k = 2.1 fm^^ (open squares), at the density of 0.34 fm~^. 

IV. MOMENTUM DEPENDENCE OF THE BHF FIELD Mi{k,uj) AT FIXED FRE- 
QUENCY 

At the density of 0.34 fm^'^, we calculate the dependence of Vi{k,u) and Wi{k,ijj) upon 
momentum k for two fixed frequencies, namely Ui = 20 MeV and uj2 = 160 MeV. The 
corresponding on-shell values of the momentum, which can be obtain from the energy- 
momentum relation u{k) = k"^ /2m + V[k^uj{k)], are approximately k{uji) = 2.15 fm^^ and 
k{co2) = 3.15 fm^^. Results for the saturation density 0.17 fm~'^ are not provided, because 
the TBF effect is not expected to play an important role at such a relatively low density. 

In Fig. [5], the upper part displays the calculated values of Vi{k, uj = 20MeV) and Vi{k, u = 
160MeV); the lower part presents the calculated values of Wi{k, u = 20MeV) and Wi{k, u = 
leOMeV). 

As we can see from the figure, for both the real and the imaginary parts of Mi{k,u), 
the open squares are very close to the corresponding filled squares in the high momentum 
region, indicating that the TBF correction is small at high momenta. However, in the low 
momentum region, the TBF has a strong effect on the shape of the fc-dependence of Vi{k, u) 
and Wi{k,u) only at the larger frequency u = 160 MeV, and it may separate the open 
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squares and the corresponding filled ones considerably. As a result, it is necessary to take 
into account the TBF effect if one wants to get more exact and reliable fc-dependence of the 
off-shell mean field Mi{k, u) felt by a nucleon with both low momentum and large frequency. 



V. SPECTRAL FUNCTION 

The spectral function S{k, u) can be calculated from Eq.(l2]), using the real and imaginary 
parts of the mass operator. Notice that in the present approximation scheme, W{k,u) = 
W2{k,LL!) for u < ep and W{k,u) = Wi{k,ijj) for u > ep. For energies u < ep, the spectral 
function S{k,ijj) is referred to as the "hole spectral function" Sh{k,u), and for energies 
CO > ep, the S{k,co) becomes the "particle spectral function" Sp{k,oj). Sh{p){k,co) measures 
the probability that a nucleon with momentum k and energy u can be removed from (added 
to) the ground state. 

In Fig. El the spectral function is plotted versus u at the density of 0.34 fm~^. The 
upper part of the figure displays the spectral distribution for momentum below the Fermi 
momentum. In the independent-particle model, states with momenta below the Fermi sur- 
face would be completely occupied so that the spectral function is identical to a 5 function 
located at the on-shell value of u. However, the two-hole configuration leads to a non vanish- 
ing imaginary part of the mass operator and consequently a finite spectral function peaked 
at the on-shell energy for momenta below kp 0, . The quasiparticle peak in the spectral 
function can be related to the shell model by the fact that when a nucleon with momentum 

k is removed from the ground state, the residual system has a large probability of having 

I I 

a well-defined excitation energy E\_i [41j. The lower part of the figure shows the spectral 
distribution for momentum above kp. 

Recently, the TBF effect on the spectral function in nuclear matter has been investigated 



46] where the Urbana TBF 
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explicitly within the in-medium T-matrix method in Ref. 
has been adopted. One may notice from the upper panel of Fig. [6] that, at momentum 
below the Fermi momentum kp, the TBF effect on the spectral distribution leads to a 
shift of the peak location to slightly higher energy and a decrease in the peak value, in 
agreement with the results of Ref. 46| within the in-medium T-matrix method using the 
Urbana TBF. It is also seen that the TBF reduces the strength of the spectral distribution 
at large negative energies. At momentum above kp, the TBF effect is mainly to shift the 
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FIG. 5: (Color online) Dependence upon k of the calculated values of Vi{k,uj) and of Wi{k^oj) for 
the two selected frequencies: w = 20 MeV and a; = 160 MeV. Open symbols correspond to the 
case with TBF contribution, while filled symbols do not. The density is fixed at 0.34 fm~^. 
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FIG. 6: (Color online) Spectral function S{k,ui) calculated from Eq.Q at the density of 0.34 fm 
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peak value to a higher energy. The TBF-induced shift of the peak location of the spectral 
distribution can be understood readily since the TBF gives an extra repulsive contribution 
to the on-shell single-particle potential and consequently increases the on-shell energy for a 
given momentum k. 

In order to test the numerical accuracy of the present work, in Fig. [7] we display the 
nucleon momentum distribution defined in Eq. (j4]) for two densities p = 0.17 and 0.34fm~^. 
By using Eq. ([5]) we get almost the same results. Due to the nucleon-nucleon correlations, 
the s.p. hole states below kp are partly empty and the particle states above kp are partly 
occupied in the correlated ground state of nuclear matter. The depletion of the lowest 
hole state ai k = at p = 0.17fm~^ is about 16.4%, which is compatible with the previous 



predictions in Refs 



• 0,3, 



35 



37 



40 



42| . This value is also consistent with the experimental 
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FIG. 7: (Color online) Nucleon momentum distribution in symmetric nuclear matter at two den- 
sities p = 0.17 fm^^ (left panel) and p = 0.34 fm^'^ (right panel) 



result in Ref. [26|]. As discussed in Ref. j40|, inclusion of the (higher order) renormalization 
contribution M3 in the mass operator may reduce the calculated depletion from ~ 17% to 
~ 14% by using a separable AVIA interaction. It is noticed that the TBF effect is negligibly 



small at the saturation density p = 0.17fm~^, in agreement with the conclusion of Ref. 35 1 
within the correlated basis function approach by adopting the Urbana f 14 interaction plus 
an effective TBF. The TBF effect only becomes sizable at high densities well above the 
saturation density as shown in the right panel of Fig. [7] where the momentum distribution 
for p = 0.34 fm^'^ is plotted. The TBF effect is shown to enhance the depletion of the 
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hole states since the TBF may induce sufficiently strong extra short-range correlations at 
sufficiently high densities. At p = 0.34 fm~^, inclusion of the TBF may enhance the depletion 
of the zero-momentum state from ~ 17% to ~ 22%. 

VI. SUMMARY 

Within the framework of Brueckner theory extended to include a microscopic TBF, we 
have calculated the dependence of the off-shell mass operator upon the momentum k and 
upon the nucleon frequency uj. The first two terms in the hole- line expansion of the mass 
operator are taken into account. Our calculations show that the TBF effect on the values 
of Mi{k,u) for fixed momentum is only important at high densities or at frequencies far 
away from its on-shell energy at kp. However, the cj-dependence of the Pauli rearrangement 
term M2{k,u) at fixed momenta is even less affected by the TBF effect. At p = 0.34 fm~^ 
which is well above the saturation density, inclusion of the TBF may enhance the repulsion 
of V2 at a large momentum k — 2.1 im~^ above kp- We also compare the off-shell values of 
Ml at fixed momenta with its on-shell values. For fixed frequency, the /c-dependence of the 
BHF field Mi is investigated, and it is shown that it is necessary to take into account the 
TBF effect if one wants to get a more exact /c-dependence of the mean field Mi{k,u) felt 
by a nucleon with both low momentum and large frequency. The nucleon spectral function 
has been calculated. At density of p = 0.34 fm~^ well above the saturation density, the 
TBF effect shifts the peak location in the spectral function to shghtly higher energy and 
reduces slightly the peak value at low momentum below the Fermi momentum kp- The TBF 
effect on the nucleon spectral function and nucleon momentum distribution turns out to be 
neglected at the saturation density p = 0.17fm~^. It becomes sizable only at high densities 
well above the saturation density, and inclusion of the TBF leads to an enhancement of the 
depletion of the zero- momentum hole state from ~ 17% to ~ 22% at p = 0.34 fm~^. 
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